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Abstract 

We study the expectation value of the phase of the fermion determinant for Wilson lattice fermions with 
chemical potential. We use quenched SU(3) ensembles and implement a recently proposed exact dimensional 
reduction of the fermion determinant. Ensembles at several temperatures below and above the phase transition 
are studied and we analyze the role of the quark mass, the temperature, the volume and the topological sectors. 
We compare our numerical results to predictions from chiral perturbation theory. 



Final version to appear in Physics Letters B. 



1. Introductory remarks 



With the increasing amount of experimental data on the QCD phase diagram, corresponding ab-initio lattice cal- 
culations become more and more important. However, when a chemical potential is introduced, lattice simulations 
face a serious challenge, the fermion sign problem. With non-zero chemical potential /j, the fermion determinant 
det[-D(//)] is complex and cannot be directly used as a probability weight. Unless conceptually new ideas are 
developed, Monte Carlo simulations need to use various kinds of reweighting techniques. 

The severeness of the sign problem, and thus the numerical effort for a reweighting strategy, may be charac- 
terized by the expectation value 

ol2e \ _ / det[£>(/x)] \ 

- / \ detpHx)] / ' {1) 

where e 1 ® is the phase of the fermion determinant dct[Z?(ju)]. Having the extra factor of 2 in the exponent on 
the Ihs. is convenient, since the expectation value of that phase may be written as a ratio of two determinants. 

Recently the determinant phase (TJJ was addressed in several papers [I] - [6J using different analytical tools, 
such as chiral perturbation theory or random matrix models. The dependence of the sign expectation value (e ) 
on the chemical potential, the volume, the temperature, the quark/pion mass and the topological sector was 
studied. On the lattice several individual results may be spotted [7\ - [10], but a systematical analysis of the sign 
problem is still missing. 

In this paper we attempt a small step towards a more complete analysis and study the determinant phase for 
quenched ensembles in a wide range of temperatures T and chemical potential values /.i. Even in the quenched 
situation this is still a sizable task, but applying the recently proposed factorization formula [TTJ [12] we are able 
to speed up the evaluation of the determinant phase considerably. We study the dependence of (e ) on the 
parameters /i, T, the quark mass to, the volume and the topological sectors. In particular we also compare the 
behavior in the low- and high temperature phases of quenched QCD. 



2. Technicalities 

In our study we use quenched ensembles generated with the Luscher-Weisz action [13] on lattices of size N 3 xNt = 
6 3 x 4, 8 3 x 4 and 10 3 x 4. The scale was determined in [14j based on the Sommer parameter. We generated 
ensembles for a wide range of temperatures between T = 210 MeV and T = 430 MeV, with the critical 
temperature for our action determined [15j to be T c = 300(3) MeV. Our statistics is between 500 and 2000 
configurations for the smaller lattices and 100 configurations for the larger ones. An overview of our ensembles is 
given in Table[T] 

For our gauge ensembles we determined the fermion determinant for the Wilson Dirac operator D(fi) with 
chemical potential /i (the lattice spacing is set to a = 1 here), 

±4 ^ 

D{ji) XlV = S x . y - k +TVI Uv ( x ) 5x+0y . ( 2 ) 

v=±l 1 

We use the convention U— V (x) = U v {x — £>)', introduce the abbreviation v for the shift vector in the v direction, 
and k is related to the bare mass to via n = 1/(4 + to). The temporal boundary conditions for the fermions are 
anti-periodic, all other boundary conditions (gluons, spatial boundary conditions for the fermions) are periodic. 

In order to speed up the evaluation of the fermion determinant for many values of [i at the same time, we use 
the dimensional reduction formula developed in [llj. We here only very briefly sketch the idea of the construction 
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N' 4 x N T 


Pa 


a [fm] 


T [MeV] 


T/T c 


m [MeV] 


# confs. 


6 3 x 4 


7.5 


0.213 


232 


0.77 


100 


500 


8 3 x 4 


7.4 


0.234 


210 


0.70 


50, 100, 200 


2000, 2000, 1000 


8 3 x 4 


7.5 


0.213 


232 


0.77 


100 


500 


8 3 x 4 


7.6 


0.194 


254 


0.85 


100 


500 


8 3 x 4 


7.7 


0.177 


279 


0.93 


100 


500 


8 3 x 4 


7.8 


0.161 


306 


1.02 


100 


500 


8 3 x 4 


7.9 


0.148 


334 


1.11 


100 


500 


8 3 x 4 


8.0 


0.135 


364 


1.21 


100 


500 


8 3 x 4 


8.1 


0.125 


396 


1.32 


100 


500 


8 3 x 4 


8.2 


0.115 


430 


1.43 


50, 100, 200 


500, 1000, 500 


10 3 x 4 


7.5 


0.213 


232 


0.77 


100 


100 


10 3 x 4 


8.1 


0.125 


396 


1.32 


100 


100 



Table 1: Table with the parameters of our ensembles. We list the size of the lattice in lattice units, the 
gauge coupling (3c, the lattices spacing a in fm, the temperature T in MeV and as a multiple of T c , the bare 
quark mass m in MeV and the statistics. 



and refer to [111 [T2J for the technical details: Applying a decomposition of the lattice into four domains, the 
fermion determinant may be rewritten in the form 

det[D(/i)] = A detfl - H - e" NT H +1 - e^^ff-i] . (3) 

Here Aq is a factor that is essentially a product of determinants for the terms of the Dirac operator restricted to 
the four domains. This factor is real and independent of the chemical potential /i. All of the /i-dependence comes 
from the second factor which has again the form of a determinant. However, the matrices Ho, H±i live on only a 
single time slice and thus the evaluation of this second determinant is by a factor of iVf. = 4 3 = 64 cheaper than 
the evaluation of the determinant in the original formulation. The matrices Ho,H±i are made from products of 
propagators on the domains and are small enough such that they can be stored in memory. Then with ([3]) the 
determinant det[D(n)] can be evaluated for several values of [i very efficiently. Actually, since Aq cancels in ([1}, 
for the phase factor e l2e only the dimensionally reduced determinant det[l — Hq — e^ NT H + i — e~^ NT H^i] is 
needed. We typically use 16 to 26 different values for /i spaced with Afj, = 0.05 (in lattice units). 

For some of our ensembles we also evaluated the topological charge of the configurations. For that purpose 
the low lying eigenvalues of the overlap operator were computed and the topological charge was determined from 
the number of zero modes using the index theorem. 



3. Numerical results 

3.1 Qualitative dependence on mass and volume 

We begin our presentation of the numerical results with a qualitative discussion of the behavior of (e l2e ) as a 
function of the quark mass and the volume. In the Ihs. plot of Fig. Q] we show (e l2e ) for the 8 3 x 4 ensemble at 
gauge coupling fie = 7.4 where we have three values of the bare quark mass, m = 50, 100 and 200 MeV. For all 
three values we observe a Gaussian-type of distribution. The width of the distribution depends on the mass with 
the largest mass giving the widest distribution, in other words small mass increases the sign problem. 
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P G = 7.4, 8 3 x4 P G = 7.5, m = 100MeV 




H [MeV] n [MeV] 

Figure 1: Phase of the determinant as a function of the chemical potential for low temperature. In the lhs. 
plot we compare the dependence on the quark mass in a fixed volume, while in the rhs. plot the mass is held 
fixed and we vary the lattice size. The symbols are connected to guide the eye. 

Similarly in the rhs. plot of Fig. [1] we compare the three volumes which we have at /3g = 7.5 for m = 100 
MeV. Here it is obvious, that the distribution is wider for small volume, i.e., increasing the volume makes the sign 
problem more severe. 

We remark that although we here only show plots for ensembles in the low temperature phase, we confirmed 
that the same behavior is also found in the high temperature phase, i.e., small mass and large volume increase 
the sign problem. This subsection is only meant to qualitatively demonstrate the behavior of (e l2e ) as a function 
of the quark mass and the volume. However, at least in the low temperature phase one can go beyond that and 
compare the numerical results to quantitative analytical predictions. This will be done in the next section. 

In this section we conclude with a check of the reliability of our numerical results. Since the phase of the 
determinant is a difficult to measure quantity - after all it is exponentially suppressed for large [i - such a 
consistency check is important. In principle the phase e i2e is a complex number. However, due to the symmetry 
under time reflections (e l2e ) is real. Thus a check for the quality of the Monte Carlo result for (e l2e ) is to inspect 
the imaginary part of that quantity and to control whether it is compatible with zero. In Fig. [2] we show the 
imaginary part of the phase expectation value for two ensembles in the low- and the high temperature phases for 
three different quark masses. In all cases we find that the imaginary part is zero within error bars as expected. 

3.2 Comparison to chiral perturbation theory 

After the first round of a more qualitative assessment of (e l2e ) in the last subsection, we now focus on the low 
temperature regime where analytical results are available. Since for our ensembles we have m^L larger than 5 for 
all our ensemblefl we can compare to results obtained in the p-regime of chiral perturbation theory [S [2], while 

1 The lattice extent L is given by L = aN, with a and N listed in Table [T] The pion masses were obtained from the fits 
to chiral perturbation theory as discussed in this subsection. 
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P G = 7.4, 8 3 x4 



P G = 8.2, 8 3 x4 




Figure 2: Imaginary part of the determinant phase as a function of the chemical potential. We show results 
for a low- and a high temperature ensemble on 8 3 x 4 at three different quark masses. The symbols are 
connected to guide the eye. 

we cannot expect agreement with the random matrix theory results in the microscopic regime. 

The basis for our comparison are the chiral perturbation theory results for (e l2e ) discussed in [2]. In particular 
the result for the expectation value of the determinant phase in the quenched case reads 

(e ae ) = e-soM + Mo) ( ) = _JZ| V — K 2 (n m v /3) cosh(n 2p0) . (4) 

n— 1 

Here V is the 4-volume, (3 = 1/aNx is the inverse temperature, and K2 denotes the modified Bessel function. 
The expression (j4|) is expected to describe the phase expectation value for /1 < m w /2. 

In Fig.[3]we show again the results for the low temperature ensembles already used in Fig. [I] but now plot them 
as a function of fj, in lattice units (which is more convenient for our fits). In the Ihs. plot we performed individual 
one parameter (the pion mass m OT ) fits for all three quark masses. The fit parameters for the three ensembles with 
bare quark masses m = 50, 100, 200 MeV, are am, = 0.812(25), 0.861(25) and 0.938(34), which corresponds to 
physical pion masses of m n = 682(21), 724(21) and 789(29) MeVd. Only data points with /i < m n /2, where (0 
is expected to hold, were taken into account in the fit. It is obvious that the corresponding curves in the Ihs. plot 
of Fig. [3] describe the data very well. Actually the values of x^/d.o.f. are rather small (below 0.1) throughout. 

The data in the rhs. plot of Fig. fallow for a much more stringent test of the chiral perturbation theory result 
(f4|. The three ensembles used there differ only by their volume, while the gauge coupling (3g and the quark mass 
m were held fixed. Thus one expects that the three ensembles have the same pion mass (neglecting possible finite 
size effects). Since the volumes V are known for the three ensembles, one can attempt a single one parameter 
(the pion mass m^) fit for all three volumes simultaneously. The result of this fit is shown by the full curves 
in the rhs. plot of Fig. [3] and obviously represents the data pretty well. The outcome for the fit parameter is 
am, = 0.729(52) which corresponds to a pion mass of — 675(49) MeV (note that here the gauge coupling 

2 We remark at this point, that the quark masses m we quote are bare quark masses for a non-chiral lattice Dirac operator 
(Wilson operator), and one cannot expect m\ oc m here. 
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P G = 7.4, 8 3 x4 P G = 7.5, m = 100MeV 




a p. a n 

Figure 3: Phase of the determinant plotted as a function of /j, in lattice units for low temperature. In the 
lhs. plot we compare the dependence on the quark mass in a fixed volume, while in the rhs. plot the mass 
is held fixed and we vary the lattice size. The symbols are our numerical data, and the full curves represent 
the fits to the chiral perturbation theory results. 

is different from the one in the lhs. plot). We find \ 2 /d.o.f '. — 0.485 which demonstrates that (j4]) describes the 
data for the different volumes very well. It is also interesting to note that the result for the pion mass from the 
fit (m^ = 675 MeV) is rather close to what one would estimate from the quenched spectroscopy calculation [16] 
done for the Luscher-Weisz action (m, = 660 MeV). 

3.3 The role of the topological charge 

An interesting effect has been discussed in |6J: In a random matrix calculation in the microscopic regime it 
was found that the phase expectation value (e ) should depend on the topological charge Q of the gauge 
configurations. The distribution of (e 126 ) becomes wider as \Q\ is increased, in other words, the sign problem is 
milder for higher charge sectors. The question which we address here is whether the topological effect is a specific 
feature of the microscopic regime or plays a role in general. 

As already remarked, we determine the topological charge via the number of zero modes of the overlap Dirac 
operator using the index theorem. The overlap operator has a free parameter s which may be used to tune the 
locality [17]. This parameter, however, also influences the number of zero modes as it shifts the center of the 
overlap projection. We compare two different values s = 0.5 and s — 0.8 which give slightly different results for 
the topological charge. The agreement of the topological charges becomes better as 0g is increased. 

In Fig. |4] we show (e l2e ) as a function of the chemical potential in lattice units for our (3g = 7.4, 8 3 x 4 
ensemble, separating the gauge ensembles with respect to the topological charge Q. The lhs. plot is for s = 0.5, 
the rhs. for s — 0.8. The plots clearly indicate that there is no topological effect for (e ) in our data. The curves 
are on top of each other within error bars, and the central values do not show the expected monotony (upwards 
trend with increasing \Q\). The same absence of the topological effect was also found for the other values of 
/3q- One may speculate about the reason for the absence of the topological effect. It could well be that it is 
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P G = 7.4, 8 3 x4, m = 100 MeV , s = 0.5 P Q = 7.4 , 8 3 x 4 , m = 1 00 MeV , s = 0.8 




a p. a u. 

Figure 4: Phase of the determinant as a function of /i in lattice units for our j3a — 7.4, 8 3 x 4 ensemble. 
In this plot we separate the data with respect to the topological charge Q of the gauge configurations. We 
compare s = 0.5 (lhs. plot) and s = 0.8 (rhs.). 

indeed seen only in the microscopic regime or only for a chirally symmetric Dirac operator such as the overlap, 
but of course also the analysis could be improved by using larger and finer lattices, where the topological charge 
is better defined. Since the cost for evaluating (e ) rises tremendously with the volume, such a study must be 
left for the future. 

3.4 Increasing the temperature 

We now discuss the behavior of (e ) as one increases the temperature. For the quenched theory the system 
undergoes a phase transition at T c ~ 300 MeV, where the center symmetry is broken spontaneously. This 
breaking is, e.g., signaled by the Polyakov loop which vanishes below T c and has a finite value above T c . In 
the high temperature phase the system spontaneously selects one out of three possible center sectors, which are 
distinguished by the phase <\>p of the Polyakov loop, <j)p ~ 0,27r/3 or — 27r/3. The fermion determinant on the 
other hand is not invariant under center transformations and becomes very small for the two complex sectors due 
to self averaging of the canonical determinants in the fugacity expansion [12J. Only for the real center sector the 
fermion determinant remains large, which in turn is the reason that in a dynamical simulation always the real 
center sector is selected. Consequently in our analysis we here consider the real center sector. We only take into 
account those gauge configurations where the phase <pp of the Polyakov loop obeys \4>p\ < 7r/3. For the other 
sectors the fermion determinant becomes relatively small and the phase factor e l20 is ill defined. 

In Fig. [5] we give an overview of the determinant phase for the 8 3 x 4 ensembles plotted as a function of 
the chemical potential in lattice units, comparing all our values of the gauge coupling. The critical value for 
the transition is roughly given by (3g = 7.8. In Fig. [6] the same information is presented as a 3-d plot. The 
figures show clearly that qualitatively the distribution of (e l2e ) versus /i keeps the Gaussian type of shape as 
the temperature is increased. The distribution first seems to become more narrow with increasing (3q, but from 
Pg = 7.8 on widens again. However, it is important to keep in mind, that although in lattice units the volume 
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Figure 5: Phase of the determinant as a function of the chemical potential in lattice units. We present 
the results for all values of the gauge coupling on our 8 3 x 4 ensembles. For these plots we only use the 
configurations in the real center sector. 



remains fixed, the lattice spacing and thus the physical volume decrease with increasing /3g- One expects that the 
shrinking physical volume has a mildening effect for the sign problem (see Subsections 3.1 and 3.2). The main 
message of Figs. [5] and [6] thus is of more qualitative nature: At the quenched transition there is only a gradual 
change in the distribution of (e ) and we do not observe a dramatic qualitative effect for the sign problem. 



4. Discussion 

In this paper we have made a step towards a more systematical understanding of the fermion sign problem for 
lattice QCD with chemical potential. Using quenched ensembles the phase of the fermion determinant is analyzed 
for a wide range of temperatures below and above the phase transition and for several values of the chemical 
potential. We analyze the dependence of the determinant phase on the temperature, the chemical potential, the 
quark mass, as well as the topological charge and compare our results to chiral perturbation theory. 

We find that for all temperatures and values of the chemical potential the sign problem becomes harder with 
increasing volume and decreasing quark mass. In general the determinant phase has a Gaussian type of distribution 
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Figure 6: 3-D plot of the phase of the determinant as a function of the chemical potential in lattice units 
and the gauge coupling. The data are for the 8 3 x 4 ensembles restricted to the real center loop sector. 

as a function of the chemical potential for all temperatures we considered. At the deconfinement transition of the 
quenched theory we do not observe any dramatic qualitative effect for the sign problem. Concerning a possible 
dependence on the topological charge, we do not find such a topological effect in the regime we work at. However, 
here technical improvements, in particular larger and finer lattices, where the concept of topological charge is 
better defined, would be needed for a final answer on the fate of the topological effect outside the microscopic 
regime. 



Acknowledgments 

We thank Jacques Bloch, Christian Lang, Stefan Olejnik and Kim Splittorff for inspiring discussions. J. D. 
acknowledges support by the FWF DK grant number W-1203, and M. M. support from the Austrian Agency for 
International Cooperation in Education and Research OeAD. Furthermore this research was supported in part by 
the Slovak Grant Agency for Science, Project VEGA No. 2/0070/09, by ERDF OP R&D, Project CE QUTE ITMS 
NFP 262401022, and by the Center of Excellence SAS QUTE (L.L.). The numerical computations were mainly 
performed on the cluster of the Department of Complex Physical Systems (Institute of Physics, Bratislava), and 
partly at the ZID, University of Graz. 

References 

[1] K. Splittorff and J.J.M. Verbaarschot, Phys. Rev. Lett. 98 (2007) 031601 | arXiv:hep-lat7 0609076|. 



8 



K. Splittorff and J.J.M. Verbaarschot, Phys. Rev. D 77 (2008) 014514 [ar Xiv:0709. 2218 [hep-lat]]. 
M. D'Elia and F. Sanfilippo. larXiv:0904.1400l [hep-lat]. 

M. P. Lombardo, K. Splittorff and J.J.M. Verbaarschot, larXiv:0904. 21221 [hep-lat], 

K. Splittorff and J.J.M. Verbaarschot, Phys. Rev. D 75 (2007) 116003 |arXiv:hep-lat/0702011] , 

J. Bloch and T. Wettig, JHEP 0903, 100 (2009) larXiv:0812.0324l [hep-lat]]. 

D. Toussaint, Nucl. Phys. B (proc. Suppl.) 17 (1990) 248. 

A. Nakamura, Y. Sasai and T. Takaishi, AIP Conf. Proc. 756 (2005) 41 6; Y. S asai, A. Nakamura and 
T. Takaishi, Nucl. Phys. Proc. Suppl. 129 (2004) 539 |a rXiv:hep-lat/0310046] . 

P. de Forcrand and V. Laliena, Phys. Rev. D 61, 034502 (2000) | ]arXiv:hep-lat/990 7004|. 



S. Ejiri, Phys. Rev. D 73 (2006) 054502 | arXiv:hep-lat/ 0506023|. 

J. Danzer and C. Gattringer, Phys. Rev. D 78 (2008) 114506 [arXiv:0809.2736 [hep-lat]]. 



C. Gattringer and L. Liptak, arXiv:0906.1088 [hep-lat]. 

M. Luscher and P. Weisz, Commun. Math. Phys. 97 (1985) 59 [Erratum-ibid. 98 (1985) 433]; G. Curci, 
P. Menotti, G. Paffuti, Phys. Lett. B 130 (1983) 205 [Erratum-ibid. B 135 (1984) 516]. 

C. Gattringer, R. Hoffmann and S. Schaefer, Phys. Rev. D 65 (2002) 094503 [arXiv:hep-lat/0112024| . 

C. Gattringer , P.E.L Rakow, A. Schafer and W. Soldner, Phys. Rev. D 66 (2002) 054502 
[arXiv:hep-lat/0202009| . 

C. Gattringer et al. [BGR Collaboration], Nucl. Phys. B 677 (2004) 3 [arXiv:hep-lat/0307013| . 



P. Hernandez, K. Jansen and M. Luscher, Nucl. Phys. B 552 (1999) 363 arXiv:hep-lat/9808010 



9 



